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ABSTRACT 

The Gamma-Ray Burst 031203 at a redshift z = 0.1055 revealed a highly reddened Type Ic 
Supernova, SN 20G31w, in its afterglow light. This is the third well established case of a link 
between a long-duration GRB and a Type Ic SN. The SN light curve is obtained subtracting the 
galaxy contribution and is modelled together with two spectra at near-maximum epochs. A red 
VLT grism 1501 spectrum of the SN near peak is used to extend the spectral coverage, and in 
particular to constrain the uncertain reddening, the most likely value for which is Eq^h{B — V) ~ 
1.07 ± 0.05. Accounting for reddening, SN 20031w is ~ 0.3 mag brighter than the prototypical 
GRB-SN 1998bw. Light curve models yield a ^^Ni mass of ~ O.55M0. The optimal explosion 
model is somewhat more massive (Mej ~ I3M0) and energetic {E ^ Qx 10^^ erg) than the model 
for SN 1998bw, implying a massive progenitor (40 — 50Afo). The mass at high velocity is not 
very large (I.4M0 above 30000km s""'^, but only O.IM0 above 60000km s~^), but is sufficient 
to cause the observed broad lines. The similarity of SNe20031w and 1998bw and the weakness 
of their related GRBs, GRB031203 and GRB980425, suggest that both GRBs may be normal 
events viewed slightly off-axis or a weaker but possibly more frequent type of GRB. 
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1. Introduction 

The first indication of a connection between 
long-duration Gamma-Ray Bursts (GRB) and 
Type Ic Supernovae (SNe Ic) was the discovery 
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of a bright SN Ic, SN 1998bw, in spatial and tem- 
poral coincidence with the nearby GRB 980425 
{z = 0.0085, Galama et al. 1998). The spectra 
of SN 1998bw showed broad P-Cygni lines of ele- 
ments such as Fe, Ca, and Si, indicating expansion 
velocities of ~ 0.1c. Models of the spectra and 
the light curve yielded an isotropic equivalent ex- 
pansion kinetic energy 5 x 10^^ erg (Iwamoto 
et al. 1998) and an ejected mass Mej ~ IIMq, 
suggesting that SN 1998bw was the highly en- 
ergetic explosion of a massive stellar core. The 
explosion synthesised a large mass of ^^Ni for a 
core-collapse event O.5M0). The progenitor 
of SN 1998bw was probably a very massive star 
(MzAMS ~ 4OM0), and the remnant very likely 
a black hole. SNe with broad spectral features, 
indicative of a large E, have been called "hyper- 
novae" . 

Nebular-phase spectra of SN 1998bw were dom- 
inated by a strong [O i] 6300, 6363A emission, as 
is typical of SNe Ic, but they also showed strong 
lines of [Fe ii], reflecting the high "''^Ni produc- 
tion. Another peculiarity was that the [O l] line 
was narrower than the [Fe 11] lines. This is not ex- 
pected in the spherically symmetric explosion of a 
CO core (Mazzali et al. 2001), and was interpreted 
as the result of the ejection of most ^^Ni along a 
preferential (polar) direction in an aspherical ex- 
plosion (Maeda et al. 2002), reminiscent of the col- 
lapsar model for GRBs (MacFadyen & Woosley 
1999). Doubts were raised about the association 
of SN 1998bw and GRB 980425 because of the 
unusual weakness of the GRB and the apparent 
presence of an X-ray transient inconsistent in po- 
sition with the SN (Plan et al. 2000), but the real- 
ity of this transient was disproved by observations 
with both XMM (Plan et al. 2004) and Chandra 
(Kouveliotou et al. 2004). 

The second clear case of association between 
a GRB and a SN was that of GRB 030329 and 
SN 2003dh. This was also a rather nearby GRB 
(z = 0.1687), but in this case it was a normal one. 
After extraction from the strong afterglow light, 
the SN turned out to be very similar to SN 1998bw 
(Stanek et al. 2003; Hjorth et al. 2003; Matheson et 
al. 2003; Lipkin 2004). A study of the hght curve 
and the spectra confirmed that SN 2003dh was 
also a very energetic SN Ic, with overall properties 
similar to those of SN 1998bw, although perhaps 
somewhat less extreme 3.5 x 10^^ erg, Mej~ 



8Mq, M{^^m)r^ 0.35Mq, Mzams - 55Mq] Maz- 
zah et al. 2003; Deng et al. 2005). A late-time 
spectrum of SN 2003dh also shows a narrow [O i] 
and broad [Fe 11] lines (Bersier et al. , in prep.). 

Given this evidence, it could be expected that 
the second closest GRB ever detected (GRB 031203 
at z = 0.1055) should also reveal SN signatures. 
The GRB was distinctly weaker than average, but 
stronger than 980425. Very strong reddening and 
a bright galaxy background made the detection of 
the SN very difficult. Nevertheless, a SN was even- 
tually detected photometrically (Tagliaferri et al. 
2004; Thomsen et al. 2004; Gal- Yam et al. 2004; 
Cobb et al. 2004) and spectroscopically (Malesani 
et al. 2004, hereafter MTC04). Perhaps not sur- 
prisingly, SN 20031w showed a spectrum similar 
to those of SNe 1998bw and 2003dh. 

The determination of the properties of SN 20031w 
requires reliable measurements of its light curve 
and spectra. This is made difficult by contamina- 
tion from host galaxy light. Therefore, we have 
reobserved the field after the SN faded in order to 
quantify the contribution of the host galaxy and 
subtract it out properly. 

In this paper we present newly reduced data for 
SN 20031w, derive the SN light curve and model 
it and the spectra. Using highly energetic SN Ic 
models similar to those used for SN 1998bw we 
determine the most likely reddening and derive the 
properties of SN 20031w. 

2. Late-time observations 

In order to perform the subtraction of the host 
galaxy light as accurately as possible, late-time 
photometry^^ and spectroscopy of the host galaxy 
were performed on 2004 May 22 (« 5 months af- 
ter the GRB) with the ESQ VLT-UTl (Antu) and 
F0RS2, at a time when the SN contribution to 
the total light should be negligible. The photom- 
etry is reported in Tab. 1. Data reduction and 
analysis were performed following standard pro- 
cedures, using the Eclipse and GAIA packages. 
The SN aperture photometry was calibrated using 
new field photometry of several standard fields. 
Photometric analysis shows only a slight fading 
(0.04 ± 0.02 mag in the R band) with respect to 

Based on observations performed at ESO-Paranal under 
programme ID 073.D-0255. 
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2004 March (MTC04). 

Spectroscopy was performed with grism 1501 
(6000 - 10000 A). The net exposure time was 45 
minutes, and the seeing was 0.6". The data were 
reduced using standard IRAF routines. This new 
"galaxy only" spectrum was used to remove the 
host contribution from a red 1501 grism spectrum 
acquired on 2003 Dec 20, which was not presented 
in MTC04 because of uncertainties in host galaxy 
subtraction. Flux calibration was achieved observ- 
ing two standard stars, LTT3864 and LTT3218. 
Consistent results were obtained for the two cases. 
In order to ensure relative calibration of the spec- 
tra of December 2003 and May 2004 we computed 
the flux of the galaxy emission lines in the differ- 
ent observations and introduced a small correction 
to make them match exactly. The entire set of 
VLT spectroscopic observations of SN 20031w is 
reported in Tab. 2. 

3. The bolometric light curve of SN 20031w 

In order to derive the properties of SN 20031w 
through modelling, it is necessary to construct the 
bolometric light curve from the observations. This 
requires careful subtraction of the host galaxy con- 
tribution. Photometry was in fact dominated by 
the host light, so that the shape and brightness of 
the reduced SN LC are sensitive to the adopted 
host brightness. Reliable host photometry can 
only be obtained when the SN has faded suffi- 
ciently. Our new data were obtained sufficiently 
late that we can assume this is the case. 

We adopted as host galaxy magnitudes the R 
and / measurements obtained on 2004 May 22 (see 
Table 1). These fluxes are both lower than those 
measured in March in these bands (MTC04), but 
the difference is at most 2a, indicating that the SN 
contributed by no more than 3% to the integrated 
flux in these bands already in 2004 March (as a 
comparison, SN 1998bw became 1 mag fainter over 
the same time interval, Patat et al. 2001). 

As for the host F-band magnitude, we had to 
base our estimate on the last available V-band 
measurement, which was taken only ~ 3 months 
after the SN explosion {V = 20.54 ± 0.05, 2004 
March 2) and is probably still affected by the SN. 
We therefore considered that the SN+host mag- 
nitude declined by 0.2 mag over the 2 months be- 
tween late Dec 2003 and 2004 Mar 2. A SN Ic 



similar to SN 1998bw declines by ~ 2.2 mag over 
that time interval. This suggests that SN 20031w 
still makes a contribution of ~ 3% to the total 
light on 2004 Mar 2. Therefore we adopted a value 
F(host) ~ 20.57 ±0.05. 

Since NIR templates for SN 1998bw are not 
available, we derived a host magnitude start- 
ing from the magnitudes measured in Feb 2004 
(MTC04), and corrected them as above assuming 
that the SN makes a contribution in J, H, and 
K similar to that in the V band. Accordingly, 
we use host galaxy magnitudes J = 18.28 ± 0.05, 
H = 17.58 ± 0.04, and K = 16.69 ± 0.06. These 
values are in fair agreement with those reported 
by Prochaska et al. (2004), Gal- Yam et al. (2004), 
and Cobb et al. (2004). 

One of the major uncertainties in the cali- 
bration of the photometry of SN 20031w is the 
reddening. This is large, and may have both 
a Galactic and a host component. Schlegel, 
Finkbeiner, & Davis (1998) report for the direc- 
tion of GRB031203/SN 20031w a galactic extinc- 
tion Eg{B -V) = 1.04. Prochaska et al. (2004) 
note that this value may be highly uncertain and 
suggest a value Eg{B — V) = 0.78 as a lower limit. 

Spectrum synthesis favours a total reddening 
Eg+h{B-V) = 1.07±0.05 (see §5). This is ~ 0.1 
smaller than the value obtained by Prochaska et 
al. (2004), and is formally consistent with the 
Galactic value of Schlegel et al. (1998) and no host 
extinction. Lack of significant extinction in GRB- 
host galaxies is indeed often reported (e.g. Watson 
et al. 2006b, and references therein). However, fol- 
lowing (Prochaska et al. 2004) we adopt Eg{B — 
V) = 0.78 and Eh{B - V) = 0.25, which results 
from the difference between Eg+h{B — V) = 1.07 
and Eg{B — V) = 0.78. Given that extinction laws 
of GRB hosts are uncertain and that an SMC-type 
law may be favoured with respect to a Milky Way- 
type law (see Watson et al. 2006b), we computed 
the intrinsic extinction suffered by the SN light 
under the 2 assumptions of a MW-type (Cardelli, 
Clayton, & Mathis 1989) and an SMC-type ex- 
tinction curve (Pei 1992) in the NIR and optical 
ranges. For Eh{B — V) = 0.25 the difference in 
the amount of extinction is at most 5% in these 
wavelength ranges, which is well within our uncer- 
tainties. Therefore we adopted a MW extinction 
curve. The value of Eh{B — V) we use here is 
likely to represent the maximum possible value. 
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In order to evaluate the bolometric light curve 
of SN 20031w we then proceeded as follows. First, 
we verified that SNe 20031w and 1998bw have 
similar spectra. We constructed the dereddened 
spectral energy distribution (SED) of SN 20031w 
at every epoch corresponding to an R and/or I- 
band measurement using the available simultane- 
ous or quasi-simultaneous optical and NIR pho- 
tometry. We interpolated the dereddened SEDs of 
SN 1998bw (using E{B - V) = 0.016 from Naka- 
mura et al. (2001)) to the phases of the SEDs 
of SN 20031w, taking into account the 1.1 time 
stretch factor. The SEDs of SN 20031w are shown 
in Figure 2 together with the reconstructed SEDs 
of SN 1998bw. In order to match the two SNe in 
the optical, the SEDs of SN 1998bw were scaled 
up by a factor 1.3. This represents the difference 
in intrinsic luminosity between the two SNe. 

Since the data of SN 1998bw do not extend to 
the IR, we extrapolated linearly the red part of the 
optical SED. The agreement between the SEDs 
of the two SNe is very good, except for the K- 
band point at the first epoch. Therefore, since the 
SEDs of SN 1998bw extend to bluer wavelengths 
than those of SN 20031w, the "uvoir" bolometric 
magnitudes of SN 20031w were computed by in- 
tegrating the flux under the SN1998bw templates 
between rest frame 3600 and 22200 A. 

The integrated "uvoir" LCs of the two SNe 
are shown in Figure 3. The LC of SN 20031w is 
broader than that of SN 1998bw by a factor of 
~ 1.1, but it is very similar to it in shape and is 
more luminous by ~ 0.3 mag. This is a smaller 
diff'erence than what was quoted in MTC04 and 
Thomson et al. (2004), and it is mostly the re- 
sult of adopting a smaller reddening to SN 20031w. 
Considering the significant systematic uncertain- 
ties 0.15 mag from the reddening to both SNe, 
^ 0.10 mag from possible peculiar velocities af- 
fecting the distance to SN 1998bw) the combined 
systematic error on the relative luminosity is ~ 
0.2 mag. 

4. Modelling the SN Light Curve 

We based synthetic LC calculations on the 1- 
D SN LC synthesis code (Iwamoto et al. 2000) 
that was used to model other Type Ic SNe (e.g. 
Nakamura et al. 2001; Mazzah et al. 2002, 2003). 
The code solves the energy and momentum equa- 



tions of the radiation plus gas in the co-moving 
frame, and is accurate to first order in v/c. Elec- 
tron densities and the electron scattering opacity 
are determined from the Saha-Boltzmann equa- 
tion. We adopted the approximation proposed by 
Gomez-Gomar & Isern (1996) for the Eddington 
factors, and fitted the TOPS opacities (Magee et 
al. 1995) to find an empirical relationship between 
the Rosseland mean and the electron scattering 
opacity. The energy deposition from radioactive 
decays was calculated with a gray 7-ray transfer 
code, assuming an absorptive opacity of O.ObYg 
cm^ g"-'^ (Swartz et al. 1995). The 7-ray energy 
source function was evaluated from the latest nu- 
clear data (Firestone et al. 1999). 

To test the above assumptions as well as othcir 
modifications we made to the code, we computed 
the LC of SN 1998bw using model CO138E50, 
which fits the spectra and LC of SN 1998bw near 
the LC peak (Nakamura et al. 2001). Model 
CO138E50 has Mej= 10.2Mq and E= 5 x 10^^ 
erg. We also used the two-component model de- 
veloped by Maeda et al. (2003). This model has 
a denser core than CO138E50 to mimic the ef- 
fect of a 2-D explosion, and thus it can repro- 
duce better the late-time LC of SN 1998bw. It 
has Mej= W.IMq and E= 4.6 x 10^^ erg. We up- 
dated the SN 1998bw bolometric LC of Patat et 
al. (2001) to a distance modulus fx = 32.76, corre- 
sponding to Hq = 72kms~^ Mpc~^. Both models 
require a ■'''^Ni mass of 0.41 ± O.O5M0. The dis- 
tribution of 56 Ni is O.O9M0 below 5,000km s-\ 
O.3OM0 between 5,000 and 23,000 km s'^, and 
O.O2M0 above 23,000 km s'^ for CO138E50 and 
0.07, 0.32, and O.O2M0 in the same three zones for 
the two-component model. These values are con- 
sistent with those of Nakamura et al. (2001), once 
the difference in the adopted distance is consid- 
ered. The early portion of the LC of SN 1998bw is 
well reproduced by both models (Figure 4). How- 
ever, as expected, the two-component model fits 
the late-time LC better than CO138E50, confirm- 
ing the results of Maeda et al. (2003). 

As we showed above, the shapes of the bolo- 
metric LCs of SNe 1998bw and 20031w are very 
similar, but the latter is broader by a factor 1.1. 
It is well known that the time scale of the LC 
near the peak, r, depends on the ejected mass 
Afej, the kinetic energy E, and the opacity k as 
T oc KV2Mej3/4£-i/4 (Arnett 1982). Since the 
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difference between the two SNe is small, we based 
our calculations on model CO138E50, increasing 
its density in order to reproduce the broader LC 
of SN 20031w. Given the relation above, we in- 
creased Mej, and consequently E,hy a, factor 1.2, 
leading to Mej= 12.25M0 and £■= 6 x lO^^ gj-g. 

We took the mass and the distribution of 
ejected ^^Ni as free parameters, and assumed 
that the rest of ejecta consists of 90% O and 
10% Si by mass. This simplified composition was 
adopted because the two spectra we have avail- 
able for modelling are too early to constrain the 
composition of the bulk of the ejecta, and also 
because our model is an approximate 1-D repro- 
duction of the probably aspherical ejecta. Addi- 
tionally, we tested the two-component model of 
SN 1998bw, also scaled up in mass by a factor 1.2, 
so that the two-component model for SN 20031w 
has Mej= 12.85M0 and E= 5.5 x 10^^ gj-g. 

Modelling the light curve alone does not yield a 
unique result: simultaneous modelling of the spec- 
tra is required, and this is performed in the next 
section. An additional complication in the case of 
SN 20031w is the uncertain reddening. We used 
a total reddening Eg+h{B — V) = 1.07, derived 
from spectral models as discussed in Section 5. 

The shape of the bolometric LC of SN 20031w is 
not much affected by the choice of reddening. Our 
main model parameters, Mej and E, affect the LC 
shape and do not change for different reddenings. 
The overall brightness of course does change, but 
it can be reproduced using different amounts of 
ejected ^'^Ni. 

Figure 5 shows the bolometric LC of SN 20031w 
and the best-fitting synthetic LCs obtained scal- 
ing CO138E50 and adjusting the ^^^Ni distribution. 
The total ^^Ni mass is ~ O.54M0. As in the case of 
SN 1998bw, we adopted a distribution of ^^Ni with 
an inner zone at v < 5, 000 km s~^ and an outer 
zone between 5,000 and 23,000 km s~^, and as- 
sumed uniform ^^Ni abundance within each zone. 
The best-fitting ^^Ni distribution is O.12M0 in the 
inner zone and O.42M0 in the outer zone. Unlike 
the case of SN 1998bw, lack of early photometry 
phase makes it impossible to tell whether ^^Ni was 
mixed to the highest velocities, so we did not in- 
troduce '^^m above 23,000km s"!. 

The synthetic LC of the scaled two-component 
model is also shown in Figure 5. The '"^^'Ni mass 



in this model is also O.54M0, of which O.IOM0 
below V < 5,000 km s"^ and O.44M0 above. The 
two-component model fits better than CO138E50 
after maximum, confirming that the dense core of 
these models is efficient in absorbing 7-rays and 
influences the LC after peak. Since the 7-ray de- 
position efficiency is increased by adding the inner 
dense core, a ^^Ni mass similar to that used for the 
rescaled CO138E50 model is sufficient to sustain 
the higher integrated luminosity of the LCs of the 
two-component models. 

Prom the above studies, we cannot select be- 
tween scaled CO138E50 and the scaled two- 
component model. The synthetic LCs differ 
before maximum, but we do not know the ex- 
act shape of the pre- maximum SN 20031w LC. 
Also, this phase is very sensitive to the ^^Ni 
distribution. Moreover, the two spectra avail- 
able are too early to show evidence of a dense 
core, unlike the cases of other previously studied 
HNe (SNe 1997ef, 1998bw, and 2002ap; Mazzah, 
Iwamoto, & Nomoto 2000; Maeda et al. 2003; 
Mazzali et al. 2002), and the light curve coverage 
does not extend to sufficiently late times to apply 
the 2-component models of Maeda et al. (2003). It 
is therefore unfortunately impossible to apply the 
2D models of the spectra and the light curve as 
was possible for SN 1998bw (Maeda et al. 2002). 
We would however expect that since SN 20031w 
was similar to SN 1998bw its late time properties 
would also be best described by a highly aspherical 
explosion viewed close to the jet axis. 

5. Spectral Modelling 

Because of the brightness of the host galaxy 
it was only possible to obtain SN spectra at two 
epochs near maximum: 2003 December 20 and 30, 
corresponding to rest-frame epochs of 16 and 24 
days past explosion, respectively. Here we add the 
red 1501 grism VLT spectrum of 2003 Dec 20 to 
the two blue 300V grism spectra shown in MTC04. 
This is useful because it is less sensitive to redden- 
ing. It shows a very strong P-Cygni profile, prob- 
ably a blend of Ca 11 IR and O i 7772A, a typical 
hypernova signature (e.g. Mazzali et al. 2002). 

We computed synthetic spectra with our Mon- 
tecarlo code (Mazzali & Lucy 1993; Lucy 1999; 
Mazzali 2000) for these two epochs. Following the 
results of the light curve modelling, we used model 
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CO138E50 rcscalcd in mass by a factor / = 1.2. 
We did not use the rescaled two-component model, 
as it only differs in the innermost part, which is not 
probed by the spectra near maximum. Our models 
can be used to constrain the total reddening, and 
have some leverage on the relative Galactic and 
host reddening since the two galaxies arc separated 
hy z = 0.1055. We adopted Eg{B -V) = 0.8 as a 
lower limit, attributed the remaining reddening to 
the host, and tested different values of rest-frame 
Eh{B - V). 

In Figure 6 we show the December 20 spectrum 
and 4 synthetic spectra, computed for Eh{B — 
V) = 0.2, 0.25, 0.3, and 0.35. The main input 
parameters of the models and the properties of 
the synthetic spectra are summarised in Table 3. 
Both Galactic and host reddening were computed 
using the Cardelli, Clayton, & Mathis (1989) law. 
Obviously, the model luminosity increases with in- 
creasing reddening. In all models the composi- 
tion is dominated by oxygen (~ 70% by mass) 
and neon (~ 20%). ''^'^Ni is only ~ 2-3% by mass, 
and silicon is even less, ~ 1-2%. This composition 
is slightly different from that used for light curve 
modelling. Including both Ne and O is important 
for the spectrum synthesis, while the synthetic LC 
is not significantly affected by replacing Ne with 
O. As reddening decreases, the abundance of ^^Ni 
also decreases, while that of oxygen increases, re- 
sulting in a progressively stronger O I absorption 
near 7000A. 

All models are very similar in the blue, where 
the main absorption features are due to Fe ii, 
Co II, and Ti ii Hnes (the big trough at 4000- 
4500A), and in the V region, which is dominated 
by lines of Si ii and S ii. However, the behaviour 
in the red, where the strong Ca ii-O i P-Cygni 
line shows prominently, is rather different. The 
model with Eh{B — y) = 0.2 shows too much ab- 
sorption in this feature. This is a consequence of 
the high abundance of oxygen. In this model, in 
fact, the abundances of the heavier elements that 
give rise to the lines in the blue are lower than 
they are in models with higher reddening, while 
the abundance of oxygen is higher. The models 
with Eh{B -V)^ 0.25, 0.3, and 0.35, all repro- 
duce the Ca ii-O i absorption well, but the model 
with the highest reddening has a large luminosity, 
and consequently a large photospheric velocity, so 
that the spectral lines are too red, and is of lower 



quality. 

Models for the December 30 spectrum are 
shown in Figure 7. Because of the smaller wave- 
length coverage at this epoch, it is not easy to 
discriminate among the models, all of which look 
reasonably good. The input values are listed in 
Table 4. The abundances are similar to those used 
for the December 20 spectrum, but ^^Ni (and de- 
cay products) is reduced by a factor of ~ 2. This 
shows that more ^^Ni was produced at higher ve- 
locities, contributing to the earlier spectrum, and 
confirms the scenario of an aspherical explosion. 
It is unfortunate that the earliest spectrum is only 
at peak, otherwise we may have been able to de- 
termine more accurately the ^^Ni abundance in 
the outermost, fastest moving layers. All models 
fit the blue part of the spectrum nicely, but those 
with Eh{B -V) = 0.25 and 0.3 seem to be the 
best in the O I-Ca II region, at least if we trust 
the red end of the spectrum. As was the case for 
the December 20 models, the oxygen abundance 
is higher for models with lower reddening, causing 
the exceedingly strong O I absorption near 7200A. 
Therefore, given our assumed Galactic reddening, 
Eg{B — V) = 0.8, wc favour a value of the host 
reddening Eh{B - y) - 0.25 ± 0.05. This results 
in a total reddening of £;H+G(S-y) ~ 1.07±0.05 
when the redshift of the host is taken into account. 
Since the total reddening is the quantity that is 
best constrained by our models, our result is also 
consistent with the Schlegel et al. (1998) galac- 
tic reddening and no host extinction. With this 
value of the reddening the absolute magnitude 
of SN 20031w on 2003 Dec 20 is -19.03, which is 
consistent with our LC analysis (Table 3). 

6. Discussion 

The properties of SN 20031w are very simi- 
lar to those of the other well-observed GRB-SNe. 
Given our chosen values of the host reddening, 
the rescaled explosion model that best fits both 
the LC and the spectra of SN 20031w has an 
Er^ 6 X 10^^ erg, M^j ~ 13Mq. From the value 
of Moj, the inferred mass of the progenitor of 
SN 20031w is - 40 - 50AfQ. These values are 
similar to, although somewhat larger than those 
of SN 1998bw. SN 20031w follows the positive 
correlation between E, M('''^Ni), and progenitor 
mass that appears to hold for SNe Ic (Nomoto 



6 



et al. 2004, 2005). As in the case of SN 1998bw, 
SN 20031w is also likely to have left behind a black 
hole remnant. 

Figure 8 shows the time evolution of the photo- 
spheric velocity as determined from spectral mod- 
elling for the three well-studied GRB-SNe and for 
other broad-lined SNc Ic that did not show an ac- 
companying GRB. The three GRB-SNe have the 
highest values, SN 2003dh being somewhat higher 
than the other two. Both SNe Ic SNe 1997ef and 
2002ap have much lower photospheric velocities. 

It must be stressed that the estimated values 
of the explosion parameters of SN 20031w are nec- 
essarily approximate. As we discussed, the esti- 
mated peak brightness of the SN, and hence the 
derived ^^Ni mass, depends sensitively on the un- 
certain galaxy subtraction and the assumed red- 
dening. Moreover, the estimate of E is based 
on spherically symmetric models of the explo- 
sion. If asymmetries were taken into account, 
most likely in the form af an aspherical explosion, 
as was possible for the much better studied case 
of SN 1998bw, it is very likely that the value of 
E would be reduced. In the case of SN 1998bw 
an estimate based on a spherical model suggested 
that the true kinetic energy is smaller by about 
a factor of 5 with respect to E, yielding a true 
E ~ lO^^erg (Maeda et al. 2002). Such a study 
was made possible for SN 1998bw by the availabil- 
ity of late-time spectrophotometry, which unfortu- 
nately cannot be obtained for SN 20031w since the 
bright host galaxy outshines the SN. The relation 
between the true and isotropic equivalent energies 
depends sensitively on both the asphericity of the 
explosion and our viewing angle. If these were not 
too different for the three GRB-SNe, as is likely, 
then the true E of SN 20031w may be reduced by 
a similar factor, to ^ 1.2 x 10"''^ erg, and that of 
SN 2003dh to - 7 x lO'^^erg. Even after such a 
major reduction, the E of the three GRB-SNe re- 
mains very large compared to that of typical SNe, 
deserving them the name "hypernovae" . 

Perhaps the most remarkable aspect is that al- 
though the properties of the three GRB-SNe vary 
by at most ±30%, those of the related GRBs (7- 
and X-ray energy output) cover ~ 4 orders of 
magnitude, GRB030329 being almost as energetic 
as a normal GRB, GRB031203 a factor of 100 
lower, and GRB980425 a factor of 100 lower still. 
The durations, 7-ray fluences and relative energy 



ranges for the 3 GRBs are reported in Table 5. 
The 7-ray outputs of the 3 GRBs (Col. 7) were 
computed according to 

E=^^. (1) 
1 + z ^ ' 

In the above equation d is the distance, S-^ the 
fluence, and z is the redshift of the GRB. We as- 
sumed isotropic emission. Since the GRBs were 
detected by different instruments, their fluences 
(Col. 3) refer to different energy ranges. However, 
these ranges do not differ much, and so we did not 
reduce the 7-ray fluxes to a common rest-frame 
energy range. The properties of the GRB-SNe are 
summarized in Table 6. 

This conundrum may be interpreted in two 
ways. One is that we may be seeing an in- 
trinsically similar phenomenon under different 
viewing angles. GRB030329/SN2003dh may be 
viewed almost pole-on, while the orientation of 
GRB980425/SN1998bw may be sufhciently off- 
axis (~ 15 - 30°, Maeda et al. 2002, 2006) 
that the detected GRB is significantly weaker. 
GRB031203/SN20031w may lie somewhere in be- 
tween (Ramirez-Ruiz et al. 2005). If this is indeed 
the case, then the 7-ray properties are a very 
strong function of angle (e.g. E.^ oc 0"'*). In this 
scenario, where line-of-sight inclination is a natu- 
ral free parameter, the optical properties of the SN 
are not much influenced by this relatively small 
spread in viewing angles, since the asymmetry in 
the SN is much less pronounced than it is in the 
GRB (Maeda et al. 2002). Larger spectral dif- 
ferences require larger inclination angles and are 
best seen in the nebular phase (Mazzali et al. 2005; 
Maeda et al. 2006) 

The other possibility is that there is a disper- 
sion in the properties of the relativistic ejecta for 
SNe with otherwise very similar characteristics. 
This may well be, especially if we consider that 
the relativistic energies at play in the GRB phe- 
nomenon (~ 10'''^ - 10^^ erg) are small compared 
to the kinetic energy involved in the SN event 
(from ^ 10''''^ erg for a normal SN to ^ 10"''^ erg for 
a hypernova) and to the presumably even larger 
neutrino energies (^ 10^^ - 10^^ erg), if there is a 
bounce (Deng et al. 2005). 

However, the inferred spread of properties of 
the SNe, albeit only of the order of ~ 30%, is seen 
not only in E and Mej , but also in less orientation- 
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dependent quantities such as the mass of ^^Ni, and 
therefore it is almost certainly real. Most likely 
both intrinsic differences and orientation effects 
arc present. So, while the observed 7-ray energy of 
GRB980425 was so low that it may require an in- 
trinsically weak GRB with normal spectral proper- 
tics, the true energy may be higher than estimated 
since the GRB was viewed off-axis as suggested by 
SN nebular line profile studies. 

The case of GRB031203 is similarly uncertain: 
based on its X-ray to 7-ray flux ratio, it was re- 
ported to be an XRF (Watson et al. 2004, 2006a). 
However, Sazonov ct al. (2004) suggested that it 
is a classical GRB based on the shape of the 7- 
ray spectrum. Based on radio calorimetry and 
the absence of the signature of an off-axis after- 
glow, Sazonov et al. (2004) and Soderberg et al. 
(2004) suggest that it is an intrinsically underlumi- 
nous event. However, Ramirez-Ruiz et al. (2005) 
model the event as a standard-energy, off-axis 
GRB. Since nebular-phase spectra of SN20031w 
arc not available, it is unfortunately impossible 
to estimate the geometry of the explosion and the 
viewing angle. 

Although it is now clear that long-duration 
GRBs and energetic, broad-lined SNe Ic (hyper- 
novae) are related (Galama et al. 1998; Hjorth et 
al. 2003; Malesani et al. 2004; Podsiadlowski et al. 
2004), the physical link between the rclativistic 
event and these SNe remains uncertain. Further 
uncontroversial occurrences of GRB-SN associa- 
tion will hopefully clarify the issue. 
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Table 1 
Late-time photometry. 



Mean date (UT) 


Exposure time 


Seeing 


Filter 


Instrument 


Magnitude 


2004 May 22.02 


2x2 min 


0.6" 


R 


FORSl 


20.44±0.02 


2004 May 22.02 


2x2 min 


0.6" 


I 


FORSl 


19.40±0.04 



Table 2 
Spectroscopy of SN 2003lw. 



Mean date (UT) 


Exposure time 


Seeing 


Grism 


Instrument 


2003 Dec 20.30 


2x15 min 


0.3" 


1501 


FORS2 


2003 Dec 20.35 


2x45 min 


0.3" 


300V 


FORS2 


2003 Dec 30.30 


2x45 min 


0.5" 


300V 


FORSl 


2004 Mar 02.12 


4x30 min 


0.6" 


300V 


FORSl 


2004 May 22.00 


3x 15 min 


0.6" 


1501 


FORSl 



Table 3 

Models for the December 20 spectrum 



Model 


Eh{B - V) 


logL 
[erg s~^] 


Vph 

km s~^ 


V 


Mv 


M(Bol) 


1 


0.20 


43.01 


17250 


21.88 


-19.43 


-18.82 


2 


0.25 


43.09 


18500 


21.90 


-19.58 


-19.03 


3 


0.30 


43.16 


18150 


21.88 


-19.75 


-19.20 


4 


0.35 


43.26 


19000 


21.83 


-19.96 


-19.45 



Table 4 

Models for the December 30 spectrum 



Model 


EniB - V) 


logL 
[erg s~^] 


Vph 

km s~^ 


V 


Mv 


M(Bol) 


1 


0.20 


42.85 


10000 


22.27 


-19.04 


-18.42 


2 


0.25 


42.91 


10500 


22.30 


-19.17 


-18.57 


3 


0.30 


43.00 


11500 


22.25 


-19.37 


-18.80 


4 


0.35 


43.09 


10625 


22.24 


-19.55 


-19.02 
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Table 5 
GRB Energetics 



GRB 


T 

s 


e-6 erg/cm2 


Range 
keV 


Instrument 


Reference 


Eiso 

erg 


980425 


31 


2.8±0.5 


40-700 


SAX-GRBM 


Pian et al. 2000 


4.3 X 10"'' 




23 


4.4±0.4 


24-1820 


CGRO-BATSE 


Galama et al. 1998 




030329 


63 


55.3±0.3 


2-30 


HETE-FREGATE 


Sakamoto et al. 2004 






63 


107.6±1.4 


30-400 


HETE-FREGATE 


Sakamoto et al. 2004 


6.9 X 




63 


163. Oil. 4 


2-400 


HETE-FREGATE 


Sakamoto et al. 2004 




031203 


40 


2.0±0.4 


20-200 


INTEGRAL-IBIS 


Sazonov et al. 2004 


4.9 X 10*^ 



Table 6 
Properties of GRB-SNe 



GRB/SN 


E 

erg 


M(^''Ni) 

Mq 


Mq 


MzAMS 

Me 


Reference 


GRB 980425/SN 1998bw 
GRB 030329/SN 2003dh 
GRB 031203/SN 20031w 


50 ± 5 
40 ± 10 
60 ± 10 


0.38-0.48 
0.25-0.45 
0.45-0.65 


10 ± 1 
8 ± 2 
13 ± 2 


35-45 
25-40 
40-50 


Iwamoto et al. 1998, Nakamura et al. 2001, Maeda et al. 2003 
Mazzali et al. 2003, Deng et al. 2005 
this paper 
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Days since GRB980425 



Fig. 1. - The i?-band LC of SN 20031w, galaxy-subtracted and time-stretched by a factor 0.9 (fiUed dots), 
compared to the V (triangles) and R (squares) LCs of SN 1998bw. Because of the redshift, the i?-band 
LC of SN 20031w is equivalent to rest-frame 6OOOA. The equivalent 6OOOA LC of SN 1998bw, obtained 
interpolating the V and R LCs, is shown as a dotted line. The LC of SN 20031w has been shifted in flux by 
an arbitrary amount (a factor 900) to match the corresponding 6OOOA LC of SN 1998bw. The shape of the 
"time compressed" light curve of SN 20031w matches well that of SN 1998bw. 
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Fig. 2. — Broad-band spectra (SEDs) of SN20031w corrected for Galactic and host reddening. These were 
computed using Hq = 72 km Mpc~^ (Spergel et al. 2003) and a standard cosmology (J^a = 0.72, Q,m = 
0.28), which yields a distance to SN 20031w of 490 Mpc (i.e. a distance modulus fi = 38.44). Only for the fluxes 
corresponding to the first epoch (top spectrum) the y-axis units are physical (erg s^^ A^^). The subsequent 
spectra, each represented by a different symbol, have been spaced apart by factors of 10, for clarity. The 
solid curves show the UBVRI SEDs of SN 1998bw, dereddened for E{B - V) = 0.016 (Nakamura et al. 
2001), reconstructed at the phases corresponding to the epochs of the SN 20031w spectra and multiplied by 
a constant factor of 1.3 to match the SN20031w spectra. The dashed curves are the IR extrapolations of the 
SN 1998bw optical spectra. 
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Fig. 3.— The "uvoir" bolometric light curves of SNe 20031w (filled dots) and 1998bw (triangles). The LC of 
SN 1998bw was computed for a distance modulus = 32.76 and E{B - V) = 0.016, that of SN 20031w for 
/X = 38.44, Eg{B -V) = 0.8, and Eh{B -V) = 0.25. The open circles show the LC of SN 20031w stretched 
by a factor of 0.9 and shifted in flux to match the LC of SN 1998bw. Error bars do not include a sistematic 
uncertainty of ±0.2 mag (Section 3) 
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Fig. 4. — Models for the "uvoir" bolometric light curve of SN 1998bw (triangles): the continuous line is 
model CO138E50, the dashed line is the 2-coniponent model. 
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Fig. 5. — Models for the "uvoir" bolometric light curve of SN 20031w (dots): the continuous line is the 
rescaled model CO138E50, the dashed line is the 2-coniponent model. 
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Fig. 6. — The VLT spectrum of SN 20031w obtained on December 20, 2003 (thick line), compared to four 
synthetic spectra, computed with Ec{B — V) = 0.8 and Eh{B — V)= 0.2 (dotted hne), 0.25 (dashed hne), 
0.3 (thin continuous hne), and 0.35 (long-dashed line), respectively. 
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Fig. 7. — The VLT spectrum of SN 20031w obtained on December 30, 2003 (thick line), compared to four 
synthetic spectra, computed with Ec{B — V)= 0.8 and Eh{B — V) = 0.2 (dotted Une), 0.25 (short-dashed 
hne), 0.3 (thin continuous hue), and 0.35 (long-dashed line), respectively. 
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Fig. 8. — Evolution of the photospheric velocity in SN 20031w (stars) and in a number of broad-lined SNe Ic. 
The references for the data are: SN 1997ef: Iwamoto et al. (2000); SN 1998bw: Iwamoto et al. (1998); 
SN 2002ap: Mazzali et al. (2002); SN 2003dh: Mazzali et al. (2004). 
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